function dx = func22(t, x, wa, T) 
%   单物种有环境 湿度变化
    global w v 
    dx = zeros(3, 1); 
    H1 = 0.8; H2 = 0.4; H3 = 0.2; 
    N1 = 2000; N2 = 2000; N3 = 2000; 
    sigma21 = 0.7; sigma32 = 0.7; sigma31 = 0.56; 
    sigma12 = 0.8; sigma23 = 0.8; sigma13 = 0.64;    
    k1 = 2; k2 =2.2;  
    if(t>60) 
        wa=-0.2+0.4*sin(0.1*t);
    end
    u = water(wa); y = temp(T); 
    sigma = k1*u*(1-w)+k2*y*(1-v); 
    dx(1) = H1*x(1)*(1 - x(1)/N1 - sigma(1)); 
    dx(2) = H2*x(2)*(1 - x(2)/N2  - sigma(2)); 
    dx(3) = H3*x(3)*(1 - x(3)/N3  - sigma(3)); 
end